A Statistical Model of Magnetic Islands in a Large Current Layer 
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We develop a statistical model describing the dynamics of magnetic islands in very large current 
layers that develop in space plasma. Two parameters characterize the island distribution: the flux 
ip contained in the island and the area A it encloses. We derive an integro-differential evolution 
equation for this distribution function, based on rules that govern the small-scale generation of 
secondary islands, the rates of island growth, and island merging. Our numerical solutions of this 
equation produce island distributions relevant to the magnetosphere and corona. We also derive 
and analytically solve a differential equation for large islands that explicitly shows the role merging 
plays in island growth. 



Recent reconnection simulations have shown that guide 
field reconnection becomes bursty f]\. Electron current 
layers near the magnetic x-line lengthen and become un- 
stable to the formation of secondary magnetic islands, 
which dynamically grow in size and begin to merge. The 
existence of these magnetic islands is consistent with ob- 
servations in both the magnetotail and the magnetopause 
, and their occurrence in current layers on the sun has 
also been inferred Q . The presence of many islands could 
also account for the observed thickness of current sheets 
associated with coronal mass ejections (CMEs), which is 
far greater than that predicted by reconnection models 
employing classical or anomalous resistivity [3| . 

Magnetic island formation is of great interest because 
of its impact on reconnection rates and because both 
simulations [5] and observations 6] indicate that islands 
are correlated with highly energetic electrons. How- 
ever, the dynamics of islands in very large-scale cur- 
rent layers are not yet well-understood. On the one 
hand, even the largest particle-in-cell (PIC) simulations 
Q come nowhere near the scale sizes of these current 
layers {L ^ 4000 di in the magnetopause, L ^ 10^ di in 
the corona). On the other hand, fluid simulations that 
encompass these large scales Q fail to capture the small- 
scale dynamics of reconnection within the dissipation re- 
gion. 

In light of these limitations, we employ a novel sta- 
tistical method for treating two-dimensional magnetic is- 
lands in large current layers. We define a statistical dis- 
trubtion function which describes islands in the whole 
current layer and develop an evolution equation for that 
distribution. We then present steady-state solutions that 
show how the merger of smaller islands drives the growth 
of large islands. The model is based on rules of island 
merging based on two-dimensional particle-in-cell recon- 
nection simulations. 

We model the islands using a distribution function / 
in terms of the island's flux ip and area A. Defined in 
this manner, dN — f{ip, A)d'ipdA signifies the number of 
magnetic islands with flux in the range [ipjip + dip] and 
area in [A, A + dA]. The distribution function f{ip, A, t) 
evolves in time in the phase space of {ip, A), independent 
of a space coordinate, because it describes islands over 
the whole current sheet of length L. The two quantities 



■0 and A are sufficient to completely characterize the state 
of an island: the magnetic field strength is computed as 
B — ijj/r — Tp^Jn/A, where r = \/ Aj-K is a characteristic 
radius. 

With the distribution function defined, we now estab- 
lish rules for the behavior of magnetic islands, and for- 
mulate how these rules affect the distribution. Islands 
first form at scales between the electron de and ion di 
skin depths in the current layers near x-lines due to the 
tearing instability The evolution equation therefore 
includes a source term S'('0, A) at this scale, in our case 
of Gaussian shape. As islands convect outwards at the 
Alfven speed, a sink term models the convection of is- 
lands out of the system at the rate ca/L. 

Furthermore, the model must account for the growth 
of magnetic islands due to reconnection. PIC simula- 
tions have shown that once reconnection reaches a non- 
linear stage, the normalized reconnection rate plateaus 
at roughly e « 0.1 Given this constant reconnection 
rate and an asymptotic reconnection magnetic field of 
Bq, the magnetic flux of an island increases at the rate 
'0 — VinBx = ecaBq. Likewise, the island's characteristic 
radius r increases at a constant rate f = eca, and so the 
area increases at the rate A = 2TrrVin = 1eca\/t^A. 

We now describe the rules for island merger: the 
merger of two islands yields an island with an area A 
equal to the sum of the two initial areas and a flux 
equal to the higher of the two initial fluxes. The flux 
does not add because magnetic reconnection does not in- 
crease or decrease magnetic flux but simply changes its 
connectivity. Thus, field lines from the island with more 
flux reconnect with those from the island with less flux 
until all of the latter's flux is depleted. The simplicity 
of these rules reflect our choice of and A as the vari- 
ables defining our phase space. Fig. [T] shows results from 
a PIC simulation that demonstrates this in the simple 
case of two isolated flux bundles. Magnetic field contours 
are shown at various times during the merging process, 
with the outermost field line representing the boundary of 
the island for the purposes of computing area. Both the 
maximum fiux ip and the total area A remained nearly 
constant throughout the simulation, with variations hav- 
ing standard deviations of 2% and 7% of the mean, re- 
spectively. These merging rules reveal why the merging 
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FIG. 1: The results of a PIC simulation of islands merger, 
showing the magnetic field lines at {a.) t — Q , (b) t — 
0.3 (c) t = 4.0 and (d) t = 8.0 fi;:/, where Q.c^ is 

the ion cyclotron frequency. The smaller island initially has 
75% of the flux and 25% of the area of the larger island. 



process is energetically favorable: the dissipation of mag- 
netic energy in the reconnection process. We can write 
the magnetic energy in an island as = B^A ^ -0^. Be- 
fore the merger the energy is given by Wi ^ ipi+ip2, and 
after the merger the energy (supposing -01 < "02 without 
loss of generality) is Wf ^1^2 <Wi. 

We now develop an equation for the island distribution 
/ that parallels the derivation of the coUisional Boltz- 
mann equation. In the absence of island merging the 
number density of magnetic islands is preserved as recon- 
nection changes an island's area and flux, while merging 
(analgous to coUisional scattering in the case of the Boltz- 
mann equation) causes a jump in the local island number. 
We first write the merging velocity of two islands with 
fluxes 01, "02 and areas Ai,A2, assuming constant mass 
density p, 



V^{'lpl,Ai,xlj2,A2) 



(^2^^2)3/2(^2^^2)3/2 



(1) 

where the dependence on Ai and A2 is implicit in ri and 
r2. For islands larger than dt, this formula resembles 
the hybrid Alfven velocity for asymmetric reconnection, 
= BiB2/4:TTp [iOjIill- For smaller islands, down to the 
electron skin depth de, « ■0iV'2c'i /47rpr^r| in agree- 
ment with the dispersion relation v^h oc fc of the whistler 
dynamics that dominate this regime. Lastly, w — > as 
r ^ 0. 
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Now we consider the number of islands AA^ 
formed by merging 
from an island of flux -0 that has merged with another is- 
land with flux 0' < 0;. Likewise, an island with area 
A must come from an island with A' < A that has 
merged with an island of area A — A' . The probability of 
those two islands merging within a time At is given by 
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' vAt 
dA'fii;,A') / d^' — f{i^',A-A')A^AA. 
Jo ^ 

(2) 

A similar analysis yields the number of islands lost 
through merging AA'^|gj^]^(0, A). Such an island is 
lost if it merges with any island of finite flux ip' and 
area A'. Again, the probability of merger depends on 
v{'tp,A,ip\A')At/L, and 

nOG poo \i 

^^Isnk = - / d"^' di^' — fi^', A')f{i^, A)A^AA . 
Jo Jo ^ 

(3) 

The rate of change of / due to merging is then calcu- 
lated as (AA^Igj.^ + AN\^^^)/Ail}AAAt. Combining this 
expression with the change in / due to reconnection, we 
get our evolution equation: 
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(4) 



On the left-hand side, the time-derivative of / results 
from the change in flux and area A due to reconnec- 
tion, which act to increase the island size and therefore 
conserve total island number. On the right-hand side, we 
have a source S{ip, A), a sink representing island convec- 
tion out of the system, and the merging terms. Consis- 
tent with the merging rules, the merging terms in Eq. (j4]) 
preserve total area, which can be shown by multiplying 
by A and integrating over A and ip. 

We first consider the case where the merging terms 
are negligible. This is valid if there are few islands in 
the layer {S{ijj,A) is small). For simplicity, we change 
variables from the area A to a characteristic radius r, 
with A — 7rr^ (although the islands need not be circular) 
and the new distribution function. 



F(0;,r) 
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/(^,A)— =27rr/(0-,^r2). 



(5) 



The evolution equation, without merging, becomes 



dF d ,■ , d , , 



^(V7,r)-^F. 



(6) 
the 

steady-state solution Fq is given by the Green's function 
1 



With a delta function source S = — ip')S{r 



Fai^p, r, 0', r') = exp 

ecA 



eL 



X S ii^j - ij') ~ Boir - r')) (7) 

where H(r) is the Heaviside function. All islands in this 
solution have a magnetic field strength Bq, equal to the 
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FIG. 2: The steady-state distribution function Faa{tp,r) for 
(a) Sjv* = 4, (b) Sn = 40, (c) Sf/ = 400, and (d) 5^* = 4000. 
The contours are logarithmic. 



ambient magnetic field. The solution int/j — r phase space 
is a decaying exponential along the line ip — Bqt, start- 
ing at tp = ip' and r — r' . The corresponding solutions 
for / are shown in dashed lines in Fig. [2l A key feature 
of Eq. ([Tj) is the characteristic island size rg = sL, which 
arises from balancing the island growth rate of f = ecA 
with a system transit time of L/ca- For example, with a 
reconnection rate of e ~ 0.1 on the magnetopause where 
L ^ 30-R_B, this simple model predicts islands of size 
^ 3Re- a survey of flux transfer events along the mag- 
netopause [13, [lj| determined that typical scale sizes are 
0.5 Re X 2 Re, for which r « I Re- Such island sizes are 
also consistent with those seen in direct observations of 
current sheets formed by CMEs [11]. 

If the merging terms are kept, the resulting integro- 
differential equation must be solved computationally. Al- 
though Eq. ^ appears to contain several parameters 
(the current sheet length L, the reconnection rate e, 
and the source term strength Sn — Jq dr dipS{ip,r) 
which gives the rate of island nucleation), in reality it 
only has one free parameter. By normalizing time to 
L/ca, island flux to sBqL, island size to eL, and the 
merging velocity to ecA, the only free parameter is the 
normalized amplitude of the source, Sj^ = Sn^L/ca- 
Although we set L — 4000 di here (a realistic value for 
the magnetopause), we can essentially model any system 
by varying Sn- 

The integration of Eq. (|3|) yields the steady state island 
distribution, Foo{ilj,r), which is shown for four values of 
Sn in Fig. O Whereas the solution without merging 
remains localized along the — B^r diagonal, the merg- 
ing terms break this symmetry: the merging process in- 
creases area but not flux so the distribution functions 
curve away from the diagonal toward larger r. One test 
of the veracity of our numerical solutions is their consis- 
tency with the moments of the evolution equation. In- 



tegrating Eq. ([4]) in ^/j and A yields an equation for the 
total number of islands N = dA d^pf{ip, A, t): 



dN 
~dt 



Sn-^N 



2L 



(8) 



where, for simplicity, we have taken the merging velocity 
to be ecA - The merging terms reduce N since two islands 
merge into one and are responsible for the term quadratic 
in A^. Eq. ([8]) can be solved analytically: 
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(9) 



(10) 



are the asymptotic number of islands in steady state and 
a time-scale to reach it. The form of N{t) predicted by 
Eq. ([9]) compares favorably to that of the numerical so- 
lution. Note from Eq. (|10p that Sn dictates the number 
of islands in the system, and so it effectively controls the 
importance of merging. For larger Sn, the merging term 
(^ N'^) in Eq. ^ dwarfs the convective loss term (~ N). 

Similarly, by multiplying Eq. ([3]) by A and integrating, 
we get a moment equation for the total area of all the 
islands in the system. At = /g°° dA d^pAf{ip, A, t): 



dAr „ CA . 

—— = ^irecATT r^T 

at L 



(11) 



where tt = dr d'tprF{ip,r,t) and the (small) 
source contribution was neglected. Since merging con- 
serves total area, the merging integrals do not contribute 
to this equation. 

Fig. [3Ka) shows the distribution of islands in radius. 



dil>Foc{'ip-,r) 



(12) 



for each Sn- The behavior of these curves for the largest 
islands can be deduced from Eq. @ for large A. In- 
tegrating over -0, assuming a constant merging velocity 
V — £ca, and expanding the integral over the merging 
terms gives an equation for f{A, t) = d^pf{ip, A, t): 



ecaAt df 



dA L ■' 



(13) 



The third term on the left side, which arises from the 
merging terms, describes how large islands grow in area 
by devouring smaller islands - the coefflcient ecaAt/L 
is the rate at which the total area At of all the smaller 
islands is consumed. By balancing the second and third 
terms on the left side we obtain a characteristic length 
scale above which growth via reconnection dominates 
growth via merging. Recalling that A — 27rec^r, we find 
that the transition occurs at f = At/^ttL. In steady 
state Eq. pT|) yields At = iiieLrT so f = ett- 
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FIG. 3: (a) The steady-state distribution function in r, given 
by Fca{r) as defined in Eq. (|12|l for various Sn- (b) The av- 
erage magnetic field strength B as a function of island radius 
r for various Sn- 



Equation (|T3)) admits a steady state solution, 

F^(r) = Ce-''/^^ r(r + erTY^/^-\ (14) 

where C is an arbitrary constant. For sufficiently large r, 
the exponential behavior of Eq. (fT4|) will dominate. In all 
four cases shown in Fig. [3] the distribution of islands at 
large radii agrees with the expression given in Eq. (|14p . 

Merging increases an island's area but not its flux, and 
hence leads to a decreased in-plane field strength B = 
ip/r. Since merging scales as iV^, one might surmise that 
the more islands in the system (e.g., because of larger 



Sn), the more merging takes place, and the smaller the 
magnetic field B. This hypothesis is borne out in Fig. 
[21 in that the steady-state solution for larger Sn tilts 
farther away from the ip — Bqv diagonal. This effect 
is seen in Fig. [3Kb), which shows, for various Sn, the 
average magnetic field strength as a function of island 
radius 



Bir) 



dil;F^{ip,r) 



(15) 



Comparisons of the present results with those of re- 
cent MHD simulations of large current layers |15[ are not 
possible since these simulations were limited to the early 
time behavior, where island merging plays a minimal role. 
Therefore, in future work we will benchmark this model 
using Hall MHD simulations to justify our equations for a 
current layer with many islands. Although observational 
data (for example, from THEMIS at the magnetopause) 
cannot directly measure Sn, the known distributions of 
island sizes and magnetic field strengths could, using our 
model, be used to infer it. An encouraging early result is 
a histogram of 1223 FTEs detected by THEMIS between 
2007 and 2008, which shows an exponential decrease in 
event number for island scale lengths above an i?^; p^ . 
Similarly, TRACE observations of supra-arcade down- 
flows in the solar corona during flares could also yield 
distributions of island size [iTj . 
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